Tunneling Spectroscopy on c-axis Yi-zCa^Z^CuaOr-^ Thin Films 
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Scanning tunneling spectroscopy was performed on c-axis Yi-^Ca^I^CusOr-i thin films for x= 
0, 0.05, O.f 5 and 0.20 at 4.2K. The measured spectra show main-gap, sub-gap and satellite features 
which scale similarly in energy versus Ca-doping, suggesting that they are associated with a single 
pairing energy. The data is analyzed with a multiband tunneling model which attributes the sub-gap 
features to the chain band and the satellite and main-gap features to the plane band for d x z_ 2 + 



s pairing symmetry. These results suggest that the superconductivity in Yi_ 
involves multiple bands. 
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Common to all the copper-oxide based high- 
temperature superconducting compounds are the CuC>2 
planes. The Cu02 planes are believed to be responsible 
for the wealth of phenomena observed in these materi- 
als particularly the formation of Cooper pairs with 
d x 2_ y 2 symmetry 0. YBa 2 Cu 3 7 _ 5 (YBCO) is pecu- 
liar among the cuprates in that it also has quasi-one- 
dimensional, metallic CuO chains, which may contain fi- 
nite superfluid density below the superconducting critical 
temperature (T c ) as suggested by various bulk probes 
[1) 0> Si- Since long-range order cannot be easily sus- 
tained in one-dimension, the apparent presence of super- 
fluid density in the chains suggests that the plane and 
chain bands are electronically coupled. Such coupling 
would imply that the superconductivity in YBCO is es- 
sentially multiband in nature [6l.l7l.la.l3.ll0l.llll]. and could 
conceivably affect its pairing symmetry. In fact, recent 
pair-tunneling experiments 1^, l3| have revealed a two- 
fold, d x 2_ y 2 + s pairing symmetry 3, HI in optimally- 
doped YBCO, with the d-wave node lines rotated away 
from the chain axis. 

At present it is not yet clear from these pair-tunneling 
experiments whether the observed d + s pairing symme- 
try is intrinsic to the planes [16[ , or it is an effect of cou- 
pling between the plane and chain bands @, 0, @, H, EH • 
Quasiparticle tunneling spectroscopy could help to elu- 
cidate this issue, by revealing multigap features in the 
excitation spectrum and by providing information about 
the pairing symmetry. For example, quasiparticle tun- 
neling experiments on MgB2 have revealed multiple s- 
wave gaps arising from multiband coupling 
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the case of YBCO, c-axis tunneling spectroscopy exper- 
iments have revealed multiple spectral features in addi- 
tion to the predominant ci- wave gap [19j, l20j, l2l|, |22j . Since 
multigap spectra could arise from multiband coupling, it 
would be helpful to understand the origin of these addi- 
tional spectral features. Furthermore, since both multi- 
band coupling in general [23| and pairing symmetry in 
the cuprates may vary with carrier doping [23], a de- 
tailed study of the spectral evolution with doping could 
yield important insights on the pairing in YBCO. 

In this Letter, we present a scanning tunneling spec- 



troscopy (STS) study of Ca-doped c-axis YBCO films 
at 4.2K. The measured spectra show main-gap, sub-gap 
and satellite features which scale similarly in energy as 
a function of Ca-doping, suggesting that they are associ- 
ated with a single pairing energy. The data is analyzed 
with a generic tunneling model which indicates that the 
main-gap and sub-gap features can be associated with 
the plane and chain bands respectively, while the satel- 
lite feature can be attributed to d + s splitting of the gap 
maximum in the plane band. These results suggest that 
the superconductivity in YBCO involves both the plane 
and chain bands. 

The Yi_ a; Ca a ;Ba2Cu307- ( 5 thin films measured were 
epitaxially grown by pulsed laser-ablated deposition 
(PLD) to ~50 nm thickness on {100}-oriented SrTi0 3 
substrates. The films were made from targets with nom- 
inal x — 0, 0.05, 0.15, 0.20 Ca-doping concentrations. 
Two to three films were made for each doping concen- 
tration resulting in average measured T c values of 89K, 
85K, 81K and 78K respectively. Cation substitution of 
Ca 2+ for Y 3+ in the unit cell is known to introduce ad- 
ditional holes into the Cu02 planes, thus overdoping the 
films and lowering the T c [25j . The films were transferred 
from the PLD chamber to the STS apparatus with less 
than 10 minutes exposure to air. STS was performed 
at 4.2K with either chemically-etched or mechanically- 
sheared Pt-Ir tips, by measuring the tunneling current / 
versus sample bias voltage V. The measured I- V spectra 
were numerically differentiated to obtain the conductance 
dl/dV spectra. 

Figure [1] shows representative dl/dV spectra taken at 
4.2K on an optimally-doped (x= 0) c-axis film. The spec- 
tra, staggered for clarity, were taken every 3 nm over the 
film surface, demonstrating good spatial homogeneity. A 
robust main-gap feature is observed along with a weaker 
above-gap, satellite feature indicated by arrows in the 
bottom spectrum of Fig. [TJ A weak sub-gap feature is 
also visible, becoming more pronounced with Ca-doping 
as will be discussed below. All three spectral features 
have often been seen in previous c-axis tunneling exper- 
iments on YBCO single crystals and thin films, by both 
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FIG. 1: (color online) Representative tunneling spectra taken 
on an optimally-doped c-axis YBa2Cu307~<s thin film at 4.2K. 
The spectra shown (staggered for clarity) were taken 3nm 
apart. Outer arrows indicate the satellite features. 
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FIG. 2: (color online) Representative tunneling spectra taken 
on a c-axis Yo.95Cao.o5Ba2Cu307_5 thin film at 4.2K. (a) and 
(b) were taken on different spots on the same film, showing 
some spectral variation but the same generic features and 
similar main-gap sizes. Inner and outer arrows indicate the 
sub-gap and satellite features respectively. 



planar-junction and STS geometries IjJ 2^, 21, 2^, 26[. 

Figure [5] shows two types of spectra measured on dif- 
ferent locations over a x= 0.05 Ca-doped film. Both 
spectral types exhibit the main-gap peaks, as well as the 
satellite features, while the spectrum in[2£b) also exhibit 
pronounced sub-gap features, as indicated by the arrows 
inside the main-gap. We also note that the main-gap 
peak heights are more pronounced in FigJSJa) than in 
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FIG. 3: (color online) Representative tunneling spec- 
trum taken on a c-axis Yo.85Cao.i5Ba2Cu307_{ thin film 
at 4.2K. (inset) Tunneling spectrum taken on a c-axis 
Yo.soCao.2oBa2Cu307_5 thin film at 4.2K showing one of the 
larger main-gaps measured. 



FigHfb). In general, all the Ca-doped YBCO thin films 
show spatial variation in the STS spectra over distances 
beyond ~20nm. This inhomogeneity was also observed 
in previous STS studies of similar Ca-doped YBCO thin 
films, indicating that this may be an effect of Ca substi- 
tution [27, 28 1 . Despite this spectral variation, the same 
generic features were consistently seen in our films up to 
x= 0.15. The main panel of figure [3] plots data for a 
x = 0.15 film, again showing the main-gap, sub-gap and 
satellite features, indicating their prevalence in Ca-doped 
YBCO. 

Figure [4] shows the evolution of these spectral features 
with Ca-doping. Assuming these features correspond to 
gap energies, a general trend of decreasing gap size versus 
increasing doping is observed. Quite remarkably, these 
spectral features show similar scaling in energy as a func- 
tion of doping, as indicated by the roughly constant ratios 
of either satellite to main-gap, or sub-gap to main-gap in 
the inset of Fig [J] While there are several scenarios that 
could lead to multiple peaks in the tunneling spectrum 
, this apparent energy scaling suggests that there is 
a single pairing energy responsible for all three spectral 
features. This pairing energy presumably corresponds to 
the predominant pairing interaction in the Cu02 planes. 

In order to understand how multiple spectral features 
could arise from a single pairing energy in YBCO, we 
discuss a multiband scenario in which the CuO chain 
band is coupled to the Cu02 plane band with an over- 
all d rr 2_„2+s pairing symmetry. First, if one charac- 



1 c pairing symmetry, 
terizes the order parameter in the plane by A p (fc) = 
An [cos ( k x ) — cos( k„ ) 1 /2 + A g with the ratio A s /A ~ 0.15 
[12l 1 1 31 ] , then the d-wave gap peak splits into main and 
satellite peaks Q. Second, in chain-plane coupling sce- 
narios with the pairing interaction being predominantly 
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FIG. 4: (color online) Average gap sizes, corresponding to the 
satellite (triangle), main-gap (square), and sub-gap (circle) 
features, plotted as a function of the average sample T c . Inset 
shows the average gap ratios, between satellite to main-gap 
(triangle) and sub-gap to main-gap (circle). 



in the plane, the gap function in the chain A c (fc) is 
proportional to A p (k), creating a sub-gap energy scale 
H, S E3] ■ We illustrate these points with a generic tun- 
neling model consisting of a two-dimensional plane band 
and a one-dimensional chain band without explicit ac- 
count of their coupling. In this simplified model, the 
normal-state plane dispersion is 

£ p (fc) = —2ti{cos(k x ) + cos(fcy) + 2t' cos(k x ) cos(k y ) 
+ t"(cos(2k x ) + cos(2fcj,))] - /ii (1) 

where t' — -0.1, t" — -0.25 and /ii = /j, — €\, where 
/j, is the chemical potential and t\ a constant potential 
associated with the planes. We describe the normal-state 
chain dispersion by 



i c = -2t 2 cos(k y ) - fj,2 



(2) 



where /12 = ^ — £2, and where the potential difference 
£2 — £1 determines the relative offset between the bottoms 
of the plane and chain bands. The parameters t\, ti rep- 
resent the nearest-neighbor overlap integrals in the plane 
and chain bands, respectively. We calculate the super- 
conducting dispersion for each band via the BCS relation 
^p.c(k) = £p, c (^) + Ap, c (^); where for the chain gap func- 
tion we take the simplest possible ansatz A c (/c) = A p (fc), 
which allows us to illustrate our model. In reality, the 
induced gap in the chains will have a more complicated 
structure @, however our model presented here is suffi- 
cient to capture many of the essential features predicted 
by a more thorough treatment, which calculates the cou- 
pling between the plane and chain layers explicitly [29| . 

We use the following parameter values {£i,i2>Mi!M2} 
= {120, 300, -86, -350} mcV to calculate the Fermi sur- 
faces of the plane (red) and chain (orange) bands shown 
in FigJSJa). Following the c-axis quasiparticle tunneling 
formalism described in Ref. 30], we calculate the dl/dV 



spectra for the plane (solid) and chain (dashed) bands 
using the d + s gap function introduced above and show 
the results in FigO^c). For comparison, we show the cal- 
culated spectra for pure d-wave gap function in Fig[5jd). 
Due to the anisotropic nature of both the gap function 
and the band dispersions, the Fermi surface does not gen- 
erally intersect with fc-space regions where the gap func- 
tion is maximum. Figure Efb) illustrates this scenario 
for the chain band, showing the d-wave and d + s gap 
functions projected onto the chain Fermi surface along 
k x . For the plane band, a similar scenario occurs. As 
a result of this band- structure enhanced gap anisotropy, 
the A:-space averaged gap edge measured by quasipar- 
ticle tunneling will tend to be smaller than A . This 
tendency is shown in Fig(5jc)-(f), where A is indicated 
by the dotted vertical lines. At present, our knowledge 
of the tunneling matrix elements between the tip and 
the individual plane and chain bands is incomplete [30l ]. 
thus the relative weighting of each band in a given tunnel 
junction is unknown. To illustrate a multiband situation, 
we assume a 30% plane to 70% chain band weighting and 
plot a combined spectrum in Fig(5]Je) and (f) [l7, 3 21 1. 
For the parameters used above, the d + s gap function 
produces three spectral features, while only two are pro- 
duced for the d-wave gap function. In comparing Fig[5je) 
with Fig(3j it is clear that the spectral features seen in 
our data can be attributed to a multiband scenario with 
d + s pairing symmetry. 

On the premise that these spectral features arise from 
the plane and chain bands, a few remarks should be made 
on the origin of the d + s pairing symmetry as well as 
doping evolution of both the s-wave component and the 
chain-plane coupling. First, O'Donovan and Carbotte 
[|| have shown that the self-consistently solved order pa- 
rameter pairing symmetry in a coupled plane and chain 
system is d + s. The predictions of their calculation are 
consistent with our experimental observations, and may 
in fact explain the origin of the s-wave component. 

Second, the 15% s-wave component as reported in pair- 
tunneling experiments 12, 13} is within the error bars of 
our data up to x = 0.15 Ca-doping. Beyond this doping 
level, the satellite feature is not discernible, as evidenced 
by the inset of Fig[3j showing one of the larger main- 
gaps observed onai = 0.20 film. Thus at present we 
cannot tell if the satellite feature has merged with the 
main-gap peak, or if the former has weakened to become 
indistinguishable from the spectral background. The sce- 
nario of the satellite and main-peak merging would sug- 
gest that the s-wave component has disappeared at the 
higher doping, leaving behind a purely d-wave pairing 
symmetry (Fig[5Jf)). Further studies using other exper- 
imental techniques would be necessary to elucidate this 
issue. 

Third, our data indicates that as YBCO is Ca-doped, 
the main-gap peaks are suppressed while the sub-gap 
features become more pronounced, suggesting that the 
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the evolution of the main-gap, sub-gap and satellite 
features with Ca-doping. These spectral features scale 
similarly in energy as a function of doping, suggesting 
that they are associated with a single pairing energy. Us- 
ing a multiband tunneling model, our analysis indicates 
that the sub-gap and main-gap features can be at- 
tributed to the chain and plane bands respectively, while 
the satellite feature could arise from d + s splitting of the 
gap maximum in the plane band. These results suggest 
that the superconductivity in Yi_ a; Ca :I ;Ba2Cu307_ ( 5 
involves multiple bands. 
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FIG. 5: (color online) Two-band quasiparticle tunneling 
model used to describe the gap, sub-gap and satellite fea- 
tures, (a) Fermi surface used in our model with the 2D plane 
band shown in red along with the ID chain band shown in 
orange, (b) Amplitudes of the d-wave and d + s gap functions 
projected on to the chain Fermi surface while moving along 
k x . (c) Calculated d+s spectra for the plane band (solid) 
and chain band (dashed) (d) Calculated rf-wave spectra for 
the plane band (solid) and chain band (dashed) (e) Combined 
plane and chain spectrum for d + s scenario, (f) Combined 
plane and chain spectrum for pure d-wave scenario. 



chain-plane coupling could be changing with doping. 
This observation may be consistent with recent Ra- 
man scattering data showing an increase in the chain- 
plane coupling with Ca-doping [24j. Although our sim- 
ple multiband tunneling model can explain the multiple 
spectral features we observed, a more rigorous theoretical 
treatment of the chain and plane coupling would be nec- 
essary to describe the peak-height evolution with doping 

1- 

Lastly, we note that tunneling studies of {110}- 
oriented Ca-doped YBCO thin films have observed a split 
zero-bias conductance peak (ZBCP) structure, which re- 
sembles the sub-gap features we observed on c-axis thin 
films [28]. However, the ZBCP splitting was observed 
to increase with Ca-doping 28], in contrast to how the 
sub-gap feature evolves in our c-axis data. This differ- 
ence clearly distinguishes the two junction orientations 
as probing different phenomena [221 ] . 

In summary, we have performed scanning 
tunneling spectroscopy measurements on c-axis 
Yi_ 2; Ca :I ;Ba2Cu307_5 thin films at 4.2K, to study 
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